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Abstract 

We present an explicit three-parameter class of p = jg, (—1/3 < 7 < 1), 
cosmological models admitting a two-dimensional group G2 of isometries act- 
ing on spacelike surfaces. The family is self-similar in the sense that it has a 
further homothetic vector field and it contains subfamilies of both (previously 
unknown) tilted and non-tilted Bianchi models with that equation of state. This 
is the first algebraically general class of solutions of this kind including dust in- 
homogeneous solutions. The whole class presents a universal spacelike big-bang 
singularity in the finite past. More interestingly, the case p = g/3 constitutes a 
new two-parameter inhomogeneous subfamily which can be viewed as a Bianchi 
V background with a gravitational wave travelling orthogonally to the surfaces 
of transitivity of the G2 group. This wave generates the inhomogeneity of the 
spacetime and is related to the sound waves tilting the perfect fluid. It seems 
to be the first explicit exact example of a gravitational wave travelling along a 
homogeneous background that has a realistic equation of state p = g/3. 
PACS number: 0402Jb, 04.40.Nr, 04.30.-w 



1 Introduction 

There is no need to mention that among the various methods used for the study of spa- 
tially inhomogeneous cosmological models, the research on exact solutions of Einstein 
field equations plays a crucial role. Due to the high non-linearity of the equations, the 
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exact solutions are necessary for the understanding of particular qualitative features 
that may constitute a guide in the study of general situations. Indeed, this has been 
the way in which many new kind of unexpected behaviours have been found. Of course, 
the exact solutions properties must be related and compared with the results obtained 
from other methods, such a combination usually leads to very powerful and general 
conclusions. For instance, when using the dynamical systems techniques in Cosmology 
M, some special exact solutions are shown to be asymptotic states of general classes of 
models. Exact solutions can also be compared with approximations or perturbations 
to check the validity of the involved expansions 0. Yet another example could be the 
study of the structure and appearance of singularities, which complements and sheds 
some light onto the singularity theorems and their conclusions, see || and references 
therein. 

The research on exact solutions is based on some physically reasonable restrictions 
used to simplify the Einstein equations. As an outstanding example, and with regard to 
geometrical properties, the existence of symmetries described by n-dimensional groups 
of motions G n (see [Q] and references therein) was the first assumption treated in 
a systematic Al way, see e.g. £|, |5], 0> leading to classifications of solutions as 
well as to fruitful techniques for their finding. In this sense, and with the study of 
spatially inhomogeneous cosmologies in mind, an important and particularly fruitful 
line of research during the last two decades has been the consideration of the class 
of spacetimes admitting a maximal 2-dimensional group of isometries Gi acting on 
spacelike surfaces. This line was somehow launched in |7j with a classification scheme 
for the particular Abelian case of these so-called "G% spacetimes" based solely on the 
properties and relations of the Killing vector fields. The classification was generalized 
for the non- Abelian case in ||, see || for a complete review. Among the classes defined 
in the Abelian case, the most simple subcase arises when there exists a family of surfaces 
orthogonal to the orbits of the group (it is then said that the group acts orthogonally 
transitively) and the two Killing vectors are mutually orthogonal, which implies that 
they are in fact hypersurface orthogonal, so that the metric can be cast in diagonal form 
in coordinates adapted to the Killing vectors. Focusing the attention on these diagonal 
Gi spacetimes, some additional assumptions have been made in order to simplify the 
field equations for a perfect fluid source, as for example the existence of additional 
homothetic or proper conformal symmetries, see fLOl |TTH and references therein. Let 
us recall here that there have also been general studies on orthogonally transitive G2 
cosmologies from a qualitative point of view, analyzing the autonomous system of first- 
order partial differential equations derivable from the Einstein field equations by using 
methods from the theory of dynamical systems |T^, [13], [1]] . The relations between some 
of the known explicit solutions and these theoretical studies were widely analyzed in 
IQ] and many references therein. 

Another important simplifying assumption for the perfect-fluid diagonal Gi space- 
times, which has received systematic attention, corresponds to the separability of the 
metric functions in coordinates that keep the diagonal form of the metric, called canon- 
ical coordinates. The case when these canonical coordinates that bring the metric func- 
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tions to a separate form are also adapted to the velocity vector of the fluid, that is to 



say, they are comoving coordinates too, was exhausted in |L4, |13|, [L5 , except for a very 



particular case identified in || that did not appear in |L3| accidentally. The general 
treatment of the separability in non-comoving canonical coordinates can be found in 



TEL M, where a classification for separable diagonal G 2 on S 2 perfect-fluid solutions 



was obtained depending on the number of linearly independent functions appearing 
in the metric, leading to a systematic procedure for the obtaining of solutions. The 
classification was exhausted, but not wholly solved, because once the machinery for 
the systematic derivation of solutions was established, the main effort was focused on 
finding solutions with special interest or physical relevance in order to study them in 
detail. 

Thus, for instance, an interesting solution was singled out in (named as 22BIIc) 
because of its 7-law equation of state which includes the relevant cases 7 = (dust 
models) and 7 = 1/3 (models for relativistic radiation). It is also interesting because it 
provides inhomogeneous generalizations of some Bianchi III, V and VI^ models found in 
|T7| , |TS| , see also . The particular dust solutions belonging to this family are actually 
included in one of the two classes of dust spacetimes studied in The solutions 



with non-zero 7, including those with 7 = 1/3, are new, though. In fact, the number 
of exact solutions for inhomogeneous spacetimes with ap = g/3 equation of state 
is rather scarce: as far as we know the first one appeared in the Wainwright-Goode 
family |L4| to which followed the Feinstein-Senovilla solution | 20fl , Davidson's pi| , the 



singularity- free metric of p2 |, their common generalization in the Ruiz-Senovilla class 



I3fl , and the non-diagonal p = 7^ family found by Mars and Senovilla |23|, |24 



The aim of this paper is to present the explicit family of solutions mentioned in 
the previous paragraph, as well as to perform an extensive detailed geometrical and 
physical study of its main features. The solutions constitute a three-parameter class 
of p = jg, — 1/3 < 7 < 1, cosmological models admitting a maximal G 2 acting on 
spacelike surfaces. The whole family is self-similar in the sense that it has a further 
homothetic vector field and it contains subfamilies of both tilted and non-tilted Bianchi 
models. This is the first inhomogeneous family with a 7-law equation of state having 
a free 7 which includes the 7 = case, something which may be very useful in order 
to study perturbations of the dust case. 

The structure of the paper is as follows. In Section ||] we introduce the line-element 
for the whole family in non-comoving canonical coordinates and show the ranges of the 
free parameters and the perfect-fluid variables, which are in turn expressed in terms of 
its velocity potential (subsections |2.1| and 772). The kinematical quantities of the fluid 
flow and their properties, as well as the deceleration parameter and the Weyl tensor 



are given in subsections and |2.4| . Then, in subsection |2.5| , we study the general 
symmetries of the spacetimes including the analysis of the special cases that arise, 
which include previously known exact 'non-tilted' Bianchi spacetimes together with 
some other new 'tilted' ones. Next, comoving canonical coordinates are introduced in 
subsection |2.6| , on the one hand to show that the metric is non-separable in comoving 
coordinates in general, and on the other hand to construct a half-null coordinate system 



3 



that will be used to make manifest the singularity structure and its type in subsection 



277| . The result is that the whole class presents a universal spacelike big-bang singularity 



in the finite past, which turns out to be of Kasner type j23L |26|. Similarly, the future 
asymptotic behaviour of the solutions is shown in subsection [2.8| . 

In Section ^ we present the most interesting particular subfamilies and limits of 
the general spacetime. In particular, two vacuum limits in the half-null coordinates 
are found in subsection |3.1| , providing two two-parameter families of pure gravitational 
pp-wave solutions. The p = g/3 subfamily is studied then in subsection p.2| , and it 
is given the interpretation of a Bianchi V background inhomogenized by means of a 
pure gravitational wave travelling along the direction orthogonal to the surfaces of 
transitivity of the G2 group. This gravitational wave is closely related to some acoustic 
waves which travel along and tilt the perfect fluid. This result is new in the sense 
that all previous works concerning propagation of waves in curved backgrounds were 
developed in the case of vacuum or massless minimally coupled scalar fields without 
potential (the latter includes the stiff fluid p = g case, see [p7|| for a pioneering treatment 
of the subject), see e.g. [EBL EM ISO, EHI] and references therein for a selection of main 



results, or in an anisotropic generalization of the stiff fluid in which the energy density 
equals the pressure on the direction of propagation of the waves f32j, and finally, 



in the case of electromagnetic fields, see |53]. There have also been other works on 
solutions describing exact solitonic perturbations of 7-law perfect fluid backgrounds 



34], BR 30j, but the formalism consists in translating the solutions of the Einstein field 



equations into equivalent five-dimensional massless scalar field spacetimes, and thus the 
backgrounds are severely restricted by some conditions on the matter content so that 
eventually only Friedmann-Lemaitre- Robertson- Walker (FLRW) spacetimes were used. 
Moreover, these perturbations give rise to anisotropies in the energy-momentum tensor. 
In our case, the gravitational wave is exact and travels on a spatially homogeneous but 
anisotropic background, and both the background and the resulting inhomogeneous 
spacetime satisfy the same realistic equation of state p = g/3. This is the first known 
example of such a situation. 

Finally, the dust subfamily is identified within the general classes found in 
in subsection |3l| Throughout the paper we follow the following conventions and 
notations. The metric has signature +2. denotes the Lie derivative with respect to 
the vector field v. Primes and dots will stand for derivatives with respect to x and t 
respectively. Greek indices run from to 3. We take units with 8ttG = c = 1. 



2 The models 

This section is devoted to presenting, in as compact a way as possible, the new family 
of spacetimes and their main geometrical and physical properties. 
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2.1 The line-element 



The line-element was derived using separability of the metric functions in non-comoving 
coordinates ]TJ] (more precisely, it arises as a particular case of the m = n = 2 (22BIIc) 
class as defined in H) and it can be written as follows 



ds 2 = F 2 (t, x) (-dt 2 + dx 2 ) + e^'-*) (l - e"* 



x 



with 



F(t, x) = exp 



// 



-(td 



cx 



.a + b 

and where we have defined the following constants: 
a = (4A + l)(z/ 2 -l)(z/-l), 



-fj,t 



with A being 



b= (6Xu -2X + u + l)(u 2 - 1), 

c = i/(2A + l)(z/ 2 + 1 + 2(2z/ -1){p+ 1)A), 

d = z/ 2 + 1 + 2(u + 3)Az/ 2 + 2(5z/ 2 - 4i/ + l)(u + 1)A 2 



A = I 2 - 1/4 



an auxiliary constant that will be used for the sake of simplicity. Actually, the constants 
satisfy the relation c = vid — X(a + b)) = vc that allows to cast the function F(t, x) in 
the alternative and possibly more convenient form given by 



F(t, x) = exp 



a + b 



(t — vx) 



(f % - 1 



Nevertheless, and for the sake of simplicity in some expressions, we prefer to keep the 
four constants a, b, c, d and F as given previously. 

The family of solutions has then three free parameters, I (or A), v and /x, although 
the latter simply provides the coordinate scaling. We obviously have A > —1/4, and 
we can choose \x > without loss of generality (see below). Furthermore, we must have 
a + b 7^ 0, which eventually will be equivalent to 



(5i/-3)A + f ± 0. 



2.2 The perfect fluid 

The line-element (||) is a solution of Einstein's field equations for a perfect-fluid energy- 
momentum tensor T a/3 — (g + p)u a u /3 + pg Q/3 (q + p ^ 0) whenever v is restricted by 

1 - v 2 > 0, (2) 
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which in turn implies the last condition of the previous subsection. The unit velocity 
vector field u reads then 



u 



— + v— 

dt dx 



The ranges for v and A immediately imply a > 0. One can also deduce that d > as 



follows: we have hv — Au + 1 > 0, W, and thus d > d\; 



-1/4 



(3-z/)(3-z/ 2 )/8 > 0. 



The energy density is given by 



Q 



b\ ^ 2 {\ + l) 



\a + bj (e^ - 1) ' 
and the equation of state is barotropic and obeys the gamma law 



P = 1Q, 

where 7 is given explicitly in terms of A by 

A I 2 - 1/4 
7 ~ A + 1 ~ Z 2 + 3/4' 

so that we have — 1/3 < 7 < 1. Notice that in the g + 3p = case A = — 1/4 <^> I = 
and the solutions admit a plane G3 on S2. 

From the above expressions is clear that the solutions have an initial big-bang 
singularity at t — 0, and this is why we have taken /1 > without loss of generality 
Subsection ^/7| is devoted to studying the singularity structure of the family, and in 
particular it will be shown that the big-bang singularity at t = is the only one for 
the whole family. 

The perfect-fluid region covers the entire manifold, and we have g > everywhere 
whenever a 2 — b 2 > 0, which is equivalent to 

(5i/ - 3)A + v < 0, (3) 

where we have taken into account that (a — b)/ (2(1 — z/ 2 )) = [y + 1)A + 1 > (so that 
q 7^ 0) which follows from the ranges for v and A. In fact, the previous condition (|3|) 
implies a + b = —2(1 — v 2 ) \{hv — 3)A + v\ > 0, and this ensures the fulfillment of both 
the dominant and the strong energy conditions on the whole spacetime. 

Since the fluid flow is irrotational it can be expressed as the normalized gradient of 
a scalar field, the so-called velocity potential a f26fl , that is 

u a = a )a / \J -a , 

where the commas stand for the partial derivative. Because of the gamma law equation 
of state, the corresponding energy density reads g = (—cr t a cr ' a )^ +1 ^ 2 ' 1 , and the velocity 
potential satisfies a homogeneous wave equation (non-linear whenever p 7^ g) which 
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gives the sound wave equation once it is linearized | 2~6fl . The velocity potential for the 
whole family is given by 



a 



M 2 (A + 1) 



b\ 2X+1 (2A + l)(a + 6) 



a + b 



/zca/1 — 



exp 



jlC 



(2A + l)(a + 6) 



vx) 



apart from an additive constant, so that the parameter v is nothing but the peculiar 
spatial fluid velocity —o'jb. The range given in (0) for the perfect fluid is consistent 
with this interpretation, so that given any 7 the corresponding subfamily contains all 
the possible values the ratio a' jo can achieve. Indeed, since this ratio is constant, this 
family of solutions does not follow an asymptotically velocity-dominated regime near 
the initial singularity except for the cases v — > 0, which could be seen as perturbations 
of homogeneous models (see section |2.5|). 



2.3 The kinematical quantities 

In order to compute the kinematical (and other) quantities for the fluid congruence 
defined by u, let us take the orthonormal tetrad {6 a } with a oc dx a in the above 
coordinate system {x a } = {t,x,y,z}. Of course, the vorticity of the fluid congruence 
vanishes identically. Regarding its acceleration, its non-vanishing components are 



ao = — 



a 1 



F~VA: 



v 



'1 



(4) 



so that the fluid flow does not follow geodesic trajectories at any point of the space- 
time, except for the special cases v = or A = 0, in which the acceleration vanishes 



everywhere. These special cases will be discussed later in subsection 2.5 



With respect to the expansion and the non-zero components of the shear tensor we 
have 







y/1 - v 2 (a + b) (el* - 1) 



a 2 (e* - l) + (a + 6)(A + l) 



(5) 



O"00 
0"11 



u 2 a n , 



022 + C"33 



°22 ~ C"33 



001 = —Mil, 

/x \av + b — e^aiy — 1)] 
1 - u 2 f' 2 {a + b) (e^-1) 
p-1 t* 



+ 



3(1 



x 



3V1 - v 2 (a + b) (e>* - 1) 
{[-3(d - vc) - a 2 } (e^ - l) - (a + 6)(2A - 1)} , 
1 _ 1 2^iZ [a + ub - ei*b(y - 1)] 



VT^(a + b) (e^t-1) ' 

where we have defined a 2 = d — vc + a(l — ^), which is indeed a positive constant for 
the given ranges of v and A. This can be easily deduced from its explicit expression 

a 2 = 2(1 - u 2 )(l + A) \(Au 2 -3v+ 1)A + (u 2 ) 



1 



v) 
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as every term between round brackets is strictly positive because of (0), (in fact 4z/ 2 — 
3u + 1 > 0, W), and so the less favorable case would correspond to A = —1/4, which 
gives a positive value for the term in square brackets for the valid range of v. 

Therefore, from expression (|5|) we see that the fluid congruence is expanding ev- 
erywhere, 9 > 0, starting with an unbounded value at the initial singularity t = and 
decreasing continuously from then on arriving eventually to zero as t tends to infinity 
A straightforward calculation shows also that 



al3 , , 2(\-v 2 \ 2 \\ 



-{t -> oo) 



29 2 v ' 3a 2 



- K a(l - v) - 2{d - vc)f + 3(1 - vfb 2 l 2 



so that the only case in which the solutions isotropize in the future is given by A = 1/2 
and v = (er n = 0, 022 = — 033 7^ 0), which is a 'comoving' family with an additional 
symmetry, as we will see later in subsection |2~5. 



Finally, we present the expression of the deceleration parameter q, whose general 
definition is 



so that it reads 



3 V ^ FvT^e 



u(0 = i(l + ? ), 



d — vc e M * + A a 2 e /it 
+ 



a + b (e^-1) « 2 (e^-l) + (a + 6)(A + l) 



from where it can be checked that at the singularity q(t — *■ 0) = 2. It is interesting 
to remark that, as follows from the previous expression and (^), q is independent of 
x, despite the inhomogeneity of the solutions. It should be stressed that this result 
holds in the above coordinate system, which is not comoving, and it is a simple con- 
sequence of the separability of the metric components in these coordinates. In the 



comoving coordinates adapted to the fluid flow, which will be given in subsection |2TG 
the deceleration parameter certainly depends on the corresponding spatial variable. 



2.4 The Weyl tensor 

Concerning the Weyl tensor, the non- vanishing scalars computed in the null tetrad (see 
§) k = 2~ 1 / 2 (0° - 1 ), I = 2" 1 / 2 (0° + 1 ), m = 2- 1 / 2 (0 2 + iO 3 ), are given by 



#2 
^0 



fi 2 l 



F 2 (a + b) 2 (e^ - 1) 



x 



(2(6A + c)-a)a + ((2A - 1)6 - 2d) b + 2e^(c + d - a)b 



12(a + b)F 2 {e^-iy 
fi 2 l 

2(a + b)F 2 (e"'- l) 2 



>* ((10A + 1)6 + (2A - l)a) + (a - 6)(4A + 1) 
a + 6)(2A - 1) + 2(c - d) - (2(c - d) - b - a 



s 



and therefore, the Weyl tensor does not vanish and there is no possible particularization 
to a conformally flat solution within the whole family. In particular, flat spacetime is 
not included in the family. 

The Petrov type is I at generic points for the general case with a maximal G2, and 
type D for the G3 on S2 case (/ = 0) as well as for the case v = 0, A = 2, which 
admits two additional isometries, belonging then to the class of LRS models (see next 
subsection P~5|) . 



The non-zero component of the magnetic part of the Weyl tensor with respect to u 
is given by 

= ^ ». - ^ »4, 

so that H(u) only vanishes in the type D cases: v = 0, A = 2 or I = 0. 
2.5 The symmetries. Special cases 

The line-element (|]) admits in general a 2-parameter group of isometries, generated 
by the two Killing vectors, 

dy ' ^ dz ' 

which obviously commute, so that the G2 is Abelian. Moreover, the metric admits a 
homothetic vector field given by 

as a + ay a + oz 

which satisfies 

It follows that the general family of solutions with c ^ belongs to the class of so- 



called 'tilted' inhomogeneous self-similar perfect-fluid models QTOj] : the velocity vector 
u is neither tangential nor orthogonal to the orbits of the 3-dimensional homothetic 
group H 3 generated by {£,,ff,(}i which is acting on spacelike hypersurfaces (S3). The 
algebraic structure of H3 is described by the following Bianchi types: Bianchi VI^ with 
h = — [(a — 2c)/2lb] 2 (Bianchi III is indeed included when h = —1) whenever lb ^ 0, 
and Bianchi V when lb = 0. In this last possibility, by taking into account that the 
case / = already admits a G3 on S2 group of isometries, the resulting group is an 
acting on S3 if I — 0. 

Coming to the possible particular cases with further symmetry we first find the 
already mentioned case with / = (A = —1/4), so that the equation of state has 
the form g + 3p = 0. This solution admits a G3 group of isometries acting multiply- 
transitively on spacelike plane ^-orbits, and the Petrov type is D. 

When I 0, the only possible cases with additional Killing vectors are all given 
by c = 0, in which case the line-element ([!]) admits the vector field £ given by (|6]) 
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restricted to c = 0, that is 

f = 2 tt + + 2 ^^r + -ri;( fl - 2W )*|-' 

ox a + b oy a + o cfe 

so that the if 3 becomes a G3 on S3 with the same Bianchi types as indicated above (see 
I~7| , || m and references therein). The explicit expression of c leads to two possibilities 



for c = 0: 

1. Case with v = 0. Now the velocity vector u is orthogonal to the orbits of the 
simply-transitive G 3 group, so that the resulting solutions belong to the following 
non-tilted Bianchi classes of spacetimes: 

• Bianchi V when A = 1/2. This homogeneous spacetime is a special case of 
the general Bianchi V family with p = g/3 found by Ruban in |T3 (line- 



element (9.20) in [0] with a = m 2 ). This is the only solution of the whole 
family such that the flow generated by u isotropizes in the future (see sub- 
section |2TT 



Bianchi III when A = 2. As already mentioned, this case corresponds to 
another algebraically special solution (Petrov type D), which actually admits 
a fourth Killing vector given by 

so that this case belongs to the LRS (G 4 on 53) models of class II in |36fl . 
As we said in subsection |2.3] , this is the only solution in (JI|), apart from 
the plane G3 on S2 case (1 = 0), having a purely electric Weyl tensor with 
respect to the fluid vector u, H(u) = (see fl37| ). 



• A one-parameter family of Bianchi VI^ spacetimes when A 7^ 1/2,2. This 
family is included in the class of evolving non-tilted Bianchi VI/j spacetimes 
for p = 7£> (table 9.4 in [[j]). Following the notation in [p]] (using a tilde for 
the quantities in |[|), the present families are included in the general cases 
with k = 1 4 — 37 1/2 (table 9.4 in |]), discovered by Uggla and Rosquist |3S 
up to quadratures. 

2. Case with 

l + v 2 

A 



2(1 - 2v){v + iy 

which satisfies g > and l/v'lO < 7 < 1 (for v E (—1, 1/2)). In this case, and for 
v 7^ 0, the perfect-fluid flow has a non-vanishing projection onto the G3 orbits, 
so that it constitutes a one-parameter family of exact 'tilted' Bianchi solutions. 
The free parameter can be chosen to be 7 with the restriction above and taking 
into account that the case 7 = 1/3 (and v = 0) falls onto the previous 'non-tilted' 
Bianchi V case. The Bianchi type for this 'tilted' homogeneous solutions is VI^ 
with h = — [z/ 2 (4A + 1)] _1 (including Bianchi III). 
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2.6 The comoving coordinates 



As is known 0, every perfect-fluid diagonal G2 solution can be written in comoving 
coordinates {T, X, y, z} keeping the diagonal form of the metric. By comoving coor- 
dinates we mean those such that u oc dx- The comoving coordinates will be useful 
for the study of the singularities that will be performed in subsection |2.7| . Also, this 
will prove that the use of comoving coordinates may sometimes be not well adapted to 
writing some solutions in explicit form, or even to look for them. 

The explicit change to comoving coordinates for the metric ([[]) is easily found to be 



t 



1 



(T + vX) 



1 



x 



{X + vT) 



(7) 



where the Jacobian of the change is 1. By writing fx = fi/ Vl — ^ 2 , the line-element 
becomes 

ds 2 = F 2 (T, X) (-dT 2 + dX 2 ) + eWi-OC-*) J( T , X) 

\e^ 1 -^ T -Vj(T,X)) 2l dy 2 + (e^-^-^J^X))" 21 dz 2 ) (8) 



where now 



F{T,X) =exp 



a + b 



{(d - cv)T - (c - ud)X) 



J X (T,X), 



and we have that 

J(T, X) = 1 - exp [-/2(T + z/X)] . (9) 

Then the fluid velocity vector field simply reads u = F~ l dx- As is obvious, this family 
of solutions is not separable in comoving coordinates for v 7^ 0. See for a study 
of the loss of separability when performing arbitrary coordinate changes (within the 
2-spaces orthogonal to the GVgroup orbits) which keep the diagonal form of the metric. 

As we can see, the structure of the line-element is perhaps not too cumbersome 
in comoving coordinates, but it is complicated enough so that the solutions had not 
been found until the Ansatz of separability in non-comoving coordinates was used. The 
structure shown in (|§|) may indicate some new Ansatzs providing, perhaps, generaliza- 
tions of these G2 spacetimes with p = 7^. 



2.7 The half-null coordinates. Singularity structure 

As mentioned previously, the solutions present an initial big-bang singularity at the 
spacelike hypersurface t = coming from the vanishing of the function J given in (|9]) 
for the metric Nevertheless, the form of the function F/J x may suggest that other 
singularities could be present, as for instance at x — > ±00. We are going to show 
that this is not the case, and consequently the only singularity of the solutions is the 
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reachable universal spacelike singularity at t = 0, which in the comoving coordinates 
is given by T + vX = 0. 

To that end, let us start by noticing that the coordinate ranges of (El) are in principle 
only restricted by 

T + uX>0. (10) 

Simple inspection on the expressions for the energy density and the Weyl scalars in- 
dicates that in this coordinate range the only other possible singular points would be 
those where the function F j J x vanishes, and that the singularities of both the Ricci 
and Weyl tensors coincide. At this point, it is very useful to perform the change to 
null coordinates {U, V} in the surfaces orthogonal to the orbits of the G 2 group, given 
by 

U=+=(T-X), V = ^=(T + X), (11) 
so that flip]) now becomes 

V(l + u) + U(l-u)>0. (12) 
Let us define now the following two constants 

K u = \/2-^-(l -v)(d+ c), K v = y/2-^—(l + u){d- c), 
a + o a + o 

whose fundamental property will turn out to be their positivity. Indeed, we have first 
of all that 

(1 - u)(d + c) = (1 - v 2 ) [2(3i/ - 1) 2 A 2 + (8A + \)v 2 + 1 

Since A > —1/4 and v 2 < 1 we have that z/ 2 (8A + 1) > — 1, which, together with 
a + b > 0, easily leads to k u > 0. On the other hand, the following explicit expression 

(1 + v )(d-c) = (l- v 2 ) [2A 2 (1 - v 2 ) + (4A + \)v 2 + 1 

directly shows that k v > too. The function F is now such that 

F(U, V)J-\U, V) = e^ u+KvV \ 

and the usefulness of the change is now clear since the region with F/J x — > 0, that is 
k u U + k v V — > — oo, can be reached within the range (|i~2|) only if U — > — oo or V — > — oo. 
In order to ascertain whether or not they are reachable, let us then bring them to finite 
values by making the typical coordinate change 

U = — \0g(K u u), V = — log(K v v), 

so that the range of the new null coordinates (u, v) is given, in principle, by u > 0, 
v > and the restriction coming from fll2"), which reads 

1 — y 1 + 1/ 

(k u u) (k v v) > 1. (13) 
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The line-element becomes then 



ds = -2J 2A (u,v)dudv + (K u u) d +- J(u,v) 
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(K u ii) d +cJ(u,v) dy + (K u u) d +cJ(u,v) dz > , (14) 
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where now 

J(U,V) = 1 - (/^u) -5 ^^*;)" 5 ^), 

and F/J x = (k u u)(k v v), so that the other possible singularities have been transported 
to uv — in the new coordinates. 

The point now is that the singularity at t — 0, which lies in the limit of the 
restriction given in ([13]) and has the form (for some constant A) 



u 



AV (d+c)/(d - c) , (15) 



does hide the other possible singularities at uv — in the sense that any endless past- 
directed causal curve from any point in our manifold terminates necessarily at t = 
(this is why it is called a universal big-bang singularity, see |§). In other words, uv = 
is not accessible within the physical spacetime, see figure [H 

2.7.1 Singularity type 

The three (non- vanishing) eigenvalues of the distortion tensor of the fluid congruence 
defined by u, 

being P a/3 = g a/3 + u a up the projector orthogonal to u, read as follows at the limit 
t 0: 

0i = f\, 2 = /Q + z), s = /(i-Z N 



where we have defined / = lim t ^ u/ (yl — u 2 FJ). The behaviour of the fluid congru- 
ence is described then by the scale factors /j {i : 1, 2, 3) defined by -u(log/j) = giving 
three different possibilities depending on the values of A at the limit t — > 0: 

• A < : l\ — > oo, I2, 13 — > 0. 

• A = : I2 — > 0, l\ and ^3 tend to a finite value. 

• A > : l 3 — > 00, lx, I2 — > 0. 

Therefore, the initial singularity is of cigar type whenever A 7^ and of pancake type 
for A = 0. This can be also inferred by means of the limit t — ► on the line-element in 
(|l]), which after changing (/xt) A dt = dr, redefining /x and absorbing some constants in 
y and z reads 



ds' 



= e a + bX (-dr 2 + (ut) i+^dx 2 
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Figure 1: Diagrams showing the singularity in the (u, v )-surfaces for the three possi- 
bilities (a) k v > k u , (b) k v < k u , (c) k v = k u . The whole spacetime is the product 
of these surfaces with the group orbits. As usual, null lines are at 45°. The fluid flow 
is indicated by arrowed lines in the region given by expression flTB|). Notice that the 
shown coordinates {t,x} would correspond to the Lorentzian coordinates related to 
{u,v} in the same manner as {T,X} are related to {U, V} (|ll|). The big-bang sin- 
gularity at i = 0, or equivalently ([15]), which is a spacelike hypersurface, is denoted 
by a thick curve which tends asymptotically to uv = (denoted by dashed lines). Of 
course, it is evident that uv = is hidden 'in the past' of the big-bang singularity in 
the physically meaningful spacetime (g > 0), and therefore the values u = or v = 
are unreachable. 
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In other words, the exponents pi such that k oc r Pi [|TJ (pp. 121) are {2pi = 2A/(1 + 
X),2p 2 = (1 + 20/(1 + A), 2p 3 = (1-20/(1 + A)} so that 



Pl+P2+P3= Pl+Pl+Pl = 1. 



Therefore, the singularity is of Kasner type [^5], ^6 



2.8 Future asymptotic behaviour of the solutions 

The future asymptotic behaviour is to be computed here as the limit T — > oo in 
comoving coordinates. (The same result is obtained, of course, by performing the 
limit t — > oo, but the final expressions will be eventually better related to well-known 
homogeneous spacetimes by using the comoving coordinates.) The line-element @ in 
that limit and for finite values of the spacelike coordinate X reads 



and thus the whole family tends to the perfect-fluid self-similar Bianchi I solutions 
[l], the three Pi exponents reading then p 1 = 1, p 2 = (1 — v)(a + 2bl)/(2(d — &>)), 
p 3 — (1 — v) [a — 2bl) I (2(d — cu)). The asymptotic behaviour for the only case in which 
the solutions isotropize in the future (A = 1/2, v = 0) is easily identified here as the 
flat FLRW solution with p — g/3 {p\ — p% — p 3 — 1). 



3 Interesting particular subfamilies 



In this section we present some of the solutions included in the general family which are 
of some physical interest. We have selected three types of subfamilies: the dust subfam- 
ily, which was actually considered at length in |TIJ as a particular family of a broader 
class of dust solutions; the p = g/3 family, which has a physically realistic equation of 
state for radiation-dominated epochs and, as we will see, it provides a simple example 
of how a gravitational wave can give rise to the inhomogeneization of the underlying 
perfect fluid; and the vacuum limits, which include some plane-wave spacetimes and, 
therefore, will also be relevant for the discussion about the inhomogeneization of the 
spacetime by the gravitational waves just mentioned. These three cases are treated in 
separate subsections. 



3.1 The vacuum limits 

The half-null coordinates of the previous section are very useful to find vacuum limits 
of the general solution. For the line-element as written in (0), the expression for the 
fluid velocity vector u transforms to 



u = J x (u, v) 



1 



^2(k u ii)(k v v) . 



(16) 
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and the energy density reads now 



Q — J 



-2A-1 



'a — 



;i-^ 2 (A + i) 



, 7 I i _i_ a-\-b -I , a + O 

a + °J (k u u) + 2(d +C ) i +2(d- c) 



q + 6 • 



From this it follows that v 2 = 1 provides vacuum limits in this coordinate system. 
Performing them in the form (|14|) we arrive at the following: 

1. Case with v = 1. The line-element becomes 

ds 2 = -2J 2X dudv + J [j 2l dy 2 + J~ 2l dz 2 ) 

where J = 1 + \f2fiv. This is a particular pp-wave [§] with <9 U as a null Killing 
vector. The Petrov type is N and the Weyl tensor takes the simple form 



= 4AZ/T J 



2 T-2A-2 



2. Case with v — — 1. Now the line-element reads 



ds 5 



2J 2X dudv + {n u u) a J 



K K U U) J 
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dy* 
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where we have put a = 2(1 + 4A)/(1 + 4A + 16A 2 ) and 6 = 8A/(1 + 4A + 16A 5 
and the function J becomes 



J = 1 — (k„w 



,-(5+6)/2 



This is again a pure gravitational pp-wave, and the type-N Weyl tensor has the 
only non- vanishing scalar 



AXlfi 2 



J 2X + 2 (K u u) d + b + 2 (8\ + 1) 



x 



1 - 4(16A 2 - 4A - l)(/t u n) a+b - 2(16A + 3)(k u u 



,(S+6)/2 



The previous vacuum solutions correspond to limits of the family of perfect fluid 



spacetimes of section |2.2| , where v was restricted by (0). Nevertheless, the line element 
given by (fj) contains further vacuum spacetimes arising whenever a — b = apart from 
the cases v — ±1, that is, if 

z/A + A + 1 =0, 

which together with the restriction A > —1/4 requires v > 3 or v < —1. The line- 
element for these cases reads 



ds 2 = F 2 (t,x) (-dt 2 + dx 2 ) + (ei^*) - e 



f(t+x) 



g 2 V 1 - -V — g 2 



A 2 ',],,-' + (ef(t^) _ e-t^V^d^ 2 



(17) 
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with 

F(t,x) = e-S» (ef(^)-e-f^) . 

This vacuum subfamily has in general a maximal Abelian G% acting simply transitively 
on spacelike hyper surf aces, the additional spacelike Killing vector field given by 

(5')£ +0Mh (S')l ■ 

and is of Petrov type I in general, so that it is the well-known Kasner metric 0. 




3.2 The p = g/3 subfamily. Propagation of gravitational waves 
in a homogeneous background 

This section is devoted to the radiation case 7 = 1/3 (A = 1/2). The condition (|3|) 
reads then v < 3/7. As has been shown, this family constitutes a generalization of the 
Ruban 'non-tilted' Bianchi V solution, and we are now going to prove that, in fact, 
the general family can be interpreted as the inhomogeneization, via the propagation of 
plane gravitational waves, of the mentioned Bianchi V Ruban solution. 

To that end, first of all it is necessary to choose the appropriate coordinate system in 
which the inhomogeneization will be more transparent. In this case, and despite what 
one might try at first, the natural choice is not the comoving coordinates (nor their half- 
null counterpart), but rather the original non-comoving coordinates {t,x,y,z}. The 
reason for this is that the change to comoving coordinates ([I]) depends explicitly on the 
parameter v, but this very parameter is the one defining the inhomogeneization. Thus, 
for the several different inhomogeneous metrics (which are selected by the particular 
values of v) , a different change to comoving coordinates is needed. In other words, in 
order to use a common coordinate system which is valid for the general subfamily, as 
well as for the particular Bianchi V metrics, one has to resort to the system {t, x, y, z}. 
This is an explicit example of the adequacy of using non-comoving coordinates in some 
occasions based on physical grounds. 

Now, we show how to write the line-element in a form that makes it explicit the 
homogeneous Bianchi V background and the travelling waves leading to its inhomo- 
geneization. By setting A = 1/2, the metric ([!]) can be rewritten as follows 

ds 2 = e fv+finh (-dt 2 + dx 2 ) + e 9v+p ™ h (e pv -^ Pinh dy 2 + e - w+ ^ ft »"dz 2 ) , (18) 
where 

e fv = e M _ 1; e gv - e /v-A« e pv = (l - e - "*) 3 , (19) 

correspond to the functions of the Bianchi V homogeneous solution, defined by v = 0, 
whereas 

/ 4i/(3i/+1) \, , / Au \ . 

finH - H ( (1 _ y2)(3 _ 7y) ) (t ~ PWi - I* [j^) (t-x), 
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are the functions linked to the inhomogeneities of this p = g/3 family of solutions. The 
energy density for this family takes the following expression 



g = I U 2 e -(*fv+f*H) ( 1±JL 

2 P \3-7v 

As we see, the inhomogeneity is driven by the same parameter v that 'tilts' the 
perfect fluids, that is, the peculiar spatial velocity. This inhomogeneity appears pri- 
marily in the transversal part of the metric as the function Pi n h, which depends only 
on the null coordinate t — x, so that pi n h is obviosuly a solution of the flat-space wave 
equation 

^nPiah = 0, (20) 

where O v denotes de d'Alembertian for the Minkowski (n) metric. Thus, this inho- 
mogeneity is constant at the null hypersurfaces t — x =const., or in other words, it 
propagates at the speed of light orthogonally to the surfaces of transitivity of the G2 
group of motions. We interpret this as signalling the existence of gravitational waves 
travelling in the Bianchi V background whose spacelike propagation direction is orthog- 
onal to the S2 orbits spanned by {y, z}. This interpretation seems to be in accordance 
with that given for the family of inhomogeneous non-diagonal stiff fluid solutions found 
by Wainwright and Marshman in (see also |30|), where the inhomogeneization is 



driven by an arbitrary function depending on t — x. In our case, the effect of the 
waves is a fortiori revealed also by the imprint which leave on the longitudinal part of 
the line-element, given by f in h- Actually, since f in h is functionally dependent on the 
velocity potential via 

a oc e finh ^, 

the implicit relation between the gravitational waves and the propagation of inhomo- 
geneities through the acoustic waves is manifest. The propagation of the gravitational 
and acoustic waves has thus two effects: it breaks the spatial homogeneity of the 
spacetime and at the same time tilts the velocity vector of the matter but keeping the 
perfect-fluid character of the matter content. 

Of course, all this has to be considered in a more careful way. For instance, in the 



case v = —1/3 we fall into the special cases analyzed in the final point of subsection |2~5 



so that the line-element (18) is in fact a Bianchi VI_3 spatially homogeneous solution, 



but now with p inh 7^ and f in h = ! The form of the longitudinal part is kept, but 
there appears a non-trivial wave-like inhomogeneity in the transversal part given by 
Pinh 7^ 0. It would seem that, in this case, the propagation of the waves would not give 
rise to the inhomogeneous trace left in /j n ^. However, this can be seen to arise as a 
rather exceptional situation because, by imposing f inh = but setting p inh = (3{t — x), 
one obtains a one-parameter family of Bianchi V-VI solutions for a non-perfect fluid 
matter content. The third Killing vector is given by 

C = 2 ± + (/* + (1 - + + (1 + v^) z^ z 
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The energy-momentum tensor of this family can be interpreted as a fluid with a non- 
zero energy flux in the direction of propagation of the wave-like inhomogeneity. Thus, 
in these cases, the waves seem to keep the spatially homogeneous character of the 
spacetime but breaking the perfect-fluid character of the matter. The restriction to 
the particular value j3 = — /i/4 leads to the mentioned case corresponding to v = 
—1/3, in which the matter content is a perfect fluid. In this exceptional case, the 
inhomogeneization effects of the travelling waves, together with that of the flux of 
energy and of the acoustic waves shown by the tilting of the fluid, seem altogether 
to balance in a final outcome of Pi n h which simply changes the Bianchi type of the 
solution. 

Summarizing, one could better say it is the inhomogeneity that generates the grav- 
itational waves and not the other way round. Indeed, clearly f in h ^ => pi n h ^ 0, but 
as we have just seen p in h ^ /^fmh ^ 0. This preferable point of view would state 
that the tilting (or the acoustic waves) generates the inhomogeneity in f in h which, in 
turn, is responsible for the appearance of the gravitational waves. 

The above paragraphs seriously indicate that one cannot talk about gravitational- 
wave propagation and its inhomogeneization properties in a naive manner. Therefore, 
we have tried to support the interpretation of pi n h representing gravitational waves 
propagating in the homogeneous background in two other independent ways. These 
are presented in what follows. 

The first way is the comparison with the available systematic studies for the de- 
scription of exact gravitational waves on homogeneous backgrounds. These studies 
were initially developed in P9"| |. The exact formalism presented in that reference, later 
used in f35|| , applies to perfect-fluid solutions of the Einstein field equations, as well 



as the degeneracies thereof, containing gravitational waves propagating over Bianchi 
I to VII backgrounds along an "algebraically" preferred direction. The case we are 
interested in is that of Bianchi V cosmologies which, after inhomogeneization takes 
place, become diagonal G2 perfect-fluid solutions. This case was termed as the sin- 
gle polarization (+) waves in Bianchi V, see [^, lOfl . Using the explicit form of the 



metric as given in (|T8| ) where py, gy, Pv correspond to functions of a fully general 
Bianchi V spacetime and taking Pi n h(t, x) and fi n h(t, x) to be free functions giving the 
inhomogeneous generalization, the Einstein equations for a perfect fluid can be split as 
follows: 

□„ (R e"^) = -{q - p)R e~^' e Sv+s ™\ (21) 
Dip = 0, (22) 

plus two other first order differential equations for f inh (longitudinal scale equations 
f29j). Here □ denotes the d'Alambertian, R = e^ v+Pinh describes the transverse scale 
expansion and ip = py — V^Pinh corresponds to the so called wave amplitude |29|. It 
must be noticed that, in this formalism, the function satisfying a source- free massless 
scalar field equation is py — V^Pinh, that is, the full combination of the function corre- 
sponding to the background plus the function carrying the inhomogeneity. This is so for 
the cases with a single (+) polarization, which restricts the Bianchi types to be I, III, 
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V or VI, and furthermore is in accordance with previous works in which the diagonal 
Einstein- Rosen [41] form of the line-element is taken and if) is the function appearing 
in the transverse part once the transitivity surface area element, which corresponds to 



R e~^ x in (jig) , is factorized, see |26|, |28], [31]]. Nevertheless, as we see in the pre- 
vious case, the transverse scale expansion R as defined in ]2^| is not equivalent to the 
transitivity surface area element in general. Actually, for some other Bianchi types this 
implies that the definition of "?/>" in [^, [10] differs from that in the "Einstein-Rosen 
view", and thereby the wave equation they satisfy are different. 

The function Pi n h(t, x) will actually satisfy the homogeneous wave equation Op inh = 
for those cases with g = p, that is, for stiff fluids including the vacuum and the 
minimally coupled massless scalar field solutions. This is so primarily because the 
so-called "transverse scale" equation ([H]) couples R to the matter source through the 
function g — p, so as long as this function vanishes, R e~ MX is completely independent 
of the matter content and satisfies a flat-space wave equation. One can argue then that 
R has to be in fact the "homogenous" function R = Ryif) in order to have models 



that reduce to the Bianchi spacetime when the waves reduce to zero, see |40[|. In 



other occasions, as was the case in the "Einstein- Rosen view" , the condition R ^R^ < 
everywhere was often imposed, see ]^(|. These simplifications, for the stiff fluid 



case, imply eventually that the d'Alambertian □ coincides with Ely, that is, with 
the d'Alambertian of the homogeneous background metric,[] and thus, since Op v = 
O v p v = 0, equation (|22| ) implies Op inh = 0. That is, both p v and p in h are solutions of 
the same wave equation ([22]). 

Fortunately, despite all the above, the simplification on R is not necessary and, 
moreover, for more realistic equations of state R evolves (as must be!) coupled with 
the matter. The former statement follows trivially in general because, in fllSD, the ho- 
mogeneous spacetime is recovered when p in h = 0, while the latter immediately implies 
that, in fact, □ ^ D v . In flUD they are related as follows 

□ = e -/« p v _ p. nh d t + p' mh d x ] , 

and therefore the coupling of R with the matter leads to the appearance of Pi n h, caus- 
ing the inhomogeneization of the transverse scale expansion, which in turn drives the 
coupling of the function ip with the matter through the operator □. The longitudinal 
scale equations would account then for fi n h. Notice that, as remarked before, in the 
present family (|i~8l) with fli~9D , it is fi n h, or equivalently a, that actually switches on the 
inhomogeneization. In fllgp together with (|I9| ) and because p in h is a function of t — x, 
we have that 



□/'„,/, = e finh o y p hlh - np v /\/s = v 'p= { s _ 7i/ J Pv, 

the second equality coming from (p2]). Thus, the propagation of pi n h is not that of a 
source-free massless scalar field, rather it is driven by the interaction with the back- 
ground geometry transverse part, clearly showing the non-linearity on the splitting of 

1 We are concentrating in the Bianchi V case, but similar scenarios arise in the rest of the Bianchi 
I to VII spacetimes fEfl . 



ap v /y/3 = - e ~( finh+fv) ( — 
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the background and the wave. The important thing here, however, is that the above 
equation is manifestly hyperbolic in character, as is obvious in its simpler form (|20D , 
and that the characteristic propagation speed of its solutions is the speed of light. 
Having no electromagnetic field present, the propagation of gravitational waves seems 
the best possibility. 

The second way to support our claim comes from the vacuum limits of the solutions 
presented in subsection |3.1|. Due to the non-linear interaction of the waves with the 



matter proper to General Relativity, and as we have seen in the previous paragraph 
for this particular case, it is not possible nor desirable to separate the metric into two 
linear terms representing the background and the waves respectively. However, one 
can try to annihilate completely one of the two terms, and then only the other must 
survive. Of course, by construction, if we set the wave-like part defined by pi n h to 
zero we obtain the Bianchi V, p = g/3, homogeneous background. This corresponds 
to putting v = 0. But can we also get rid of the matter, so that only the gravitational 
wave remains? A striking and beautiful answer would come from finding that the 
vacuum limits in the coordinates of flTB| ) correspond to plane gravitational waves, as 
was the case for the limits found in subsection |3.1| . Unfortunately, there are no vacuum 
limits in these coordinates. To achieve the vacuum limits one has to make use of 
other coordinate systems, but then again the limit depends on the coordinate system 
chosen |p9|| . In spite of this, the existence of coordinate systems in which the limits 



v = ±1 correspond to a gravitational pp-wave has already been shown in subsection 



3T| . At this point, the question giving unequivocal sense to the gravitational wave 
inhomogeneity interpretation for this family of p = g/3 solutions would thus be: do 
all the vacuum limits of this family with v = — 1 correspond to a plane gravitational 
wave? We do not have a rigorous answer for this question yet, but we do believe that 
the answer is positive. In this sense, we claim that getting rid of the matter, whenever 
this is possible, and within the allowed range of the parameter v, provides a pure and 
simple gravitational pp-wave, which we interpret as the remanent of the mixed case. 

Hitherto, the work on gravitational waves in cosmological backgrounds based on 
exact solutions has been mainly tackled in the cases of vacuum, scalar and electro- 
magnetic fields, and stiff fluid (and its anisotropic generalization |32| ) as sources in 



G2 on S2 spacetimes, see the reviews [2S[ ^9], |30, |31j] and references therein. As men- 



tioned in the Introduction, there are other works presenting solitonic perturbations of 
p = g/3 FLRW spacetimes later generalized to p = jg in P5| , but the energy- 
momentum tensor of the inhomogeneous spacetime turns out to be a non-perfect fluid. 
The general formalism used in these papers assumes that the backgrounds in which 
the perturbations propagate have a restricted type of energy-momentum tensors, so 
it has been applied only to FLRW spacetimes. The present family of solutions may 
constitute then the first exact solution in General Relativity for a perfect fluid with 
a realistic p = g/3 equation of state describing gravitational waves travelling on a 
spatially homogeneous background. 



21 



3.3 The dust subfamily 



The particular dust cases, defined by A = (so that requires now v < 0), belong 
to a more general family of algebraically general dust spacetimes already presented in 
fT9|| . We devote this short subsection to identifying the present dust one-parameter 



subfamily within the more general dust family appearing in |19 . 

By taking profit of the form of the line-element in the half-null coordinates {u, v, y, z} 
([14]), and by making a coordinate change to Lorentzian coordinates analogous to that 



in (|TTD, the line-element fll4]) with A = coincides exactly with expression (7) of [|T| 
for the following particular values of the constants (here, we denote with a tilde the 
quantities which appear in reference [T§[], if necessary): 



/V2. This can be seen as the free parameter of the solution. 

• C2 = 0, Ci < 0. This implies the relation k u = k v , which follows from A = 0. 

• b = v{\ + z/)/(l + z/ 2 ). By remembering that v G (—1,0), this leads to the 
restriction b G [b_, 0) for the parameter b of |19| . 

The diagram for this subfamily corresponds to the figure 2(f) of [19 with the point p 
at t = x = 0, in agreement with the particular case (c) of the present figure |l|. 
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